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Summary 

Nitrates  are  observed  in  polar  snow,  fim  and  ice  in  variable  amounts  and  concentrations. 
Hundred- year  or  longer  rime  series  of  nitrate  observations  can  be  obtained  from  ice  cores. 
It  is  important  to  identify  possible  and  probable  sources  of  nitrates  and  mechanisms  by 
which  surface  deposition  and  incorporation  into  ice  occurs.  Several  important  sources  of 
nitrate  are  solar  protons  and  auroral  electrons  producing  ionization  in  the  D-region. 
Chemical  reactions  energized  by  this  ionization  produce  nitrates.  Transport  mechanisms 
can  move  nitrates  from  the  stratosphere  into  the  troposphere  where  it  becomes  incorporated 
into  precipitation.  From  the  amounts  of  nitrates  deposited  in  polar  snow  over  long  times  it 
may  be  possible  to  determine  atmospheric  characteristics  of  importance  to  climate  and 
global  change  studies.  The  specific  goal  of  this  study  is  to  evaluate  carefully  and 
completely  the  atmospheric  nitrates  produced  by  solar  flare  protons  and  auroral  electrons 
and  to  estimate  the  amounts  deposited  in  snow. 

Statement  of  Work  and  Work  Accomplished 

There  are  several  phases  to  this  problem.  The  initial  phase,  to  refine  and  correct  the  solar 
proton  ionization  deposition  code  for  the  polar  atmosphere  is  complete  as  documented  in 
Attachment  A.  Both  the  COSPAR  and  the  MSIS  atmosphere  models  have  been  used. 
Attachment  B  describes  the  improvements  to  the  computer  code  used  by  previous 
investigators.  Attachment  C  describes  the  introduction  of  an  improved  description  of  the 
galactic  cosmic  ray  contributions  to  atmospheric  ionization.  The  development  of  a 
simplified  model  of  the  chemical  reactions  connecting  ionization  to  nitrate  production  has 
begun.  The  next  phase  of  work  will  include  the  acquisition  and  incorporation  of  a  model  of 
the  auroral  ionization  from  the  NOAA  Space  Environment  Laboratory.  Also  to  be  done  is  a 
one  and  two  dimensional  model  of  the  transport  of  nitrate.  Mr.  Vitt  is  work  this  summer 
under  the  supervision  of  Dr.  Charles  Jackman  of  the  Goddard  Space  Flight  Center  learning 
the  skills  necessary  to  run  this  model.  He  is  also  attending  the  summer  lectures  in  Earth 
Science  at  GSFC  to  acquire  further  science  background  on  atmospheric  processes. 

Publications 

At  this  time  there  are  only  internal  memos  (attachments  A,B,  and  C). 

Some  results  may  be  ready  for  publication  in  the  coming  year  and,  of  course,  the  complete 
results  will  be  published  in  Mr.  Vitt's  PhD  dissertation. 


Participating  Professionals 


Mr.  Francis  Vitt  is  being  supervised  in  his  PhD  studies  jointly  by  Profs.  Thomas 
Armstrong,  Thomas  Cravens,  Ed.  Zeller  and  Adjunct  Prof.  Gisela  Dreschoff  of  the 
Department  of  Physics  and  Astronomy.  Dr.  Charles  Jackman  of  the  Goddard  Space  Right 
Center  will  be  appointed  Adjunct  to  the  Department  of  Physics  and  Astronomy  faculty  for 
the  purpose  of  serving  as  the  external  member  of  Mr.  Vitt's  PhD  committee. 

Interactions 

Mr.  Vitt's  is  spending  6  weeks  this  summer  at  the  Goddard  Space  Right  Center  interacting 
with  science  staff  and  other  students  there. 

Discoveries  and  Applications 

The  computer  code  being  developed  will,  upon  completion,  be  available  to  the  government 
and  public  for  the  purpose  of  assisting  in  the  evaluation  of  the  global  nitrate  budget. 
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Apendix  A 


MEMO 

TO:  Charles  Jackman 

FROM:  F.  Vitt,  T.  P.  Armstr  cftgfc’Laird,  T.  E  .  Cravens,  G.  Dreschhoff,  E.  Zeller 

Referring  to  your  e-mail  of  last  year  and  our  earlier  discussions  of  the  different 
results  of  our  ionization  calculations,  we  have  examined  the  problem.  We  believe  that  we 
understand  why  the  results  differed. 

1 .  Laird  et  al.  underestimated  the  full  upper  hemisphere  flux.  Jackman  method  of 
integrating  over  incidence  angle  is  better.  Laird  method  now  includes  this. 

2.  Jackman  approximation  of  spectrum  is  fairly  coarse— overestimating  and 
underestimating  because  of  the  4  segment  fit.  Also,  our  contention  is  that  power  law  is 
much  more  suitable  for  solar  events  at  these  energies. 

The  details  are  given  below. 

Regarding  the  calculation  of  ion  pairs  produced  by  energetic  protons  incident  on  the 
upper  atmosphere  using  IMP8  data  which  has  been  done  by  Claude  Laird  and  Charlie 
Jackman,  daily  average  IMP8  data  for  day  293  of  1989  has  been  used.  Integration  of  ion 
pairs  per  cubic  cm  over  the  altitude  shows  that  Jackman's  method  tends  to  give  more 
ionization  than  Laird's  method,  by  about  a  factor  of  two.  (See  figure  1.  Here  the  circles 
represent  Laird's  results,  triangles  represent  Jackman's  results.)  To  find  the  discrepancy 
between  the  two  methods  the  following  comparisons  were  done.  Modifications  were  done 
to  the  Laird  method  to  attempt  to  obtain  results  that  match  Jackman's. 

Comparison  of  the  two  methods: 

*  Both  assume  that  the  number  of  ion  pairs  produced  equals  the  energy  deposited 
by  the  protons  in  the  atmosphere  divided  by  35  eV. 

*  Laird  uses  the  power  law  energy  spectrum  ^ =  where  is  the 

differential  flux  (no.  of  particles  *cm-2*sr'I*sec.-1*  MeV-1),  A  is  the  flux  constant,  5  is 


the  spectral  exponent.  A  and  8  are  determined  over  nine  energy  ranges  to  fit  the  IMP8 
data.  This  gives  a  piece  wise  continuous  spectrum.  (See  figure  2.) 

*  Jackman  uses  an  exponential  energy  spectrum  given  by  ^  =  ^*exP(^),  where 
F  and  Eq  are  determined  over  only  four  energy  ranges  to  fit  the  IMP8  data.  This  gives 
relatively  large  discontinuities  in  the  energy  spectrum.  (See  figure  3.) 

*  Laird's  energy  spectrum  ranges  from  0.38  MeV  to  289  MeV. 


*  Jackman's  energy  spectrum  ranges  from  0.29  MeV  to  440  MeV. 

*  Laird's  method: 

Integration  of  the  energy  spectrum  gives 


=  total  no.  of  particles  cm'2  sr1  incident  on  the  upper 
atmosphere  for  the  day, 

where  Ej  (i  =  1,2,. ..,50)  are  logarithmically  spaced  energies  ranging  from  0.38  MeV  to 
289  MeV.  This  gives  49  mono-energetic  fluxes  with  incident  energy  E?  =  VEi*Ej+1„ 


These  fluxes  are  followed  vertically  through  the  atmosphere,  which  is  divided  into 
5  km  slabs,  as  energy  deposited  in  each  slab  is  calculated,  using  a  table  look  up  method. 
Hence  calculating  energy  deposited  cm-2  s-1.  Then  Laird  multiplies  by 

UxiO^eV  wk  Ts) t0  Set  *e  units  of  no-  of  ion  cm'3- 

*  Jackman's  method: 


Jackman  assumes  an  isotropic  distribution  of  pitch  angle  and  azimuthal  symmetry. 
Jackman  calculates  no.  of  ion  pairs*cm-3*sec.-1  using  the  expression 


where  is  the  differential  flux  (no.*cm-2*sr1*sec_1*MeV-1),  w  =  35  eV  per  ion  pair, 

dE 

p(z)  is  the  atmospheric  density  at  the  altitude  z,  ^  is  the  energy  deposited  per  gm. 

Let  Ed  =  dE  be  the  energy  deposited  by  a  particle  of  energy  E  with  a  pitch  angle  a 
in  the  segment  zd,  zd  =  width  of  the  atmospheric  layer  in  gm*cm-2  centered  on  the 
altitude  z. 


>cc,z)  =  flux  of  particles  with  incident  energy  Eo  at  a  pitch  angle  a  at  the 
altitude  z  is  known,  then  qi,j(A,z)  =  ionization  rate  by  a  flux  of  particles  in  the  energy 

interval  Eq  -  to  Eo  +  and  in  the  pitch  angle  interval  a  -  to  a  +  over  the 
atmospheric  distance  zq  centered  on  the  altitude  z  can  be  calculated  with  the  expression 


qi,j(A,z)  =27t 


^<Eo,g,z) 
w  zd  seca 


sina  da  dE 


Assuming  Ed  and  ^<Eo,a,z)  change  slowly  enough  within  the  pitch  angle  interval 
dE 


we  can  write 


rB0«# 


Jr*0  2  i2 

■  ^jg<Eo,a,z)  Ed  dE  I  cos  a  sin  a 

&>-S!  /a-a 


da 


=  [sin2(a+Ap-)  -  sin2(a-  4?-)]  |  ^<Eo,a,z)  ED 


2  I  dE 


dE 


To  calculate  the  energy  deposited  Jackman  uses  the  expression  for  the  range  as 
R(E)  =  A*EB,  thus 

ED  =  E-(-^ieca +EB)b 
A 


Hence,  the  ionization  rate  at  the  altitude  z  is 


q(z)  =  p(z)  X  qij(A.z) 


=  no.  of  ion  pairs  cm~J  sec 


1-3  QPJ — 1 


Modifications  to  the  Laird  method: 

Assume  an  isotropic  distribution  of  pitch  angle  between  0  and  ^  and  azimuthal 
symmetry.  Let  a*  be  50  evenly  spaced  angles  between  0  and  ^  the  fluxes  with  pitch  angle 
between  a*  and  ock+i  have  the  pitch  angle  a®  such  that  cos  ak  =  Ycos  ak*cos  ak+i. 

Let  q;jk  =  the  no.  of  ion  pairs*cm~3  produced  by  the  flux  with  incident  energy  Ef 
in  the  j*  slab  centered  at  altitude  zj  with  pitch  angle  ak .  So, 


qi,j,k  = 


35xlO'6MeV 


.f 

L  AXj,t 


where  Eo.i.j.k  is  the  energy  deposited  by  a  proton  with  incident  energy  E° ,  pitch  angle  ak, 
in  the  j*  slab.  AXj<k  is  the  effective  thickness  of  the  j*  slab,  i.e.  AXjk  =  AZj*secak,  AZj 
is  the  thickness  of  the  slab  (5km).  Ni  (for  i  =  1,2,...49)  is  the  mono-energetic  flux  with 
incident  energy  E-* .  And  Aft  is  the  solid  angle  between  ak  and  ak+i  as  <X>  varies  from  0  to 


Assuming  Ni,  Eo.i.j.k,  and  AXj,k  do  not  change  significantly  as  pitch  angle  varies 
between  ak  and  ak+i,  we  can  write 


qi.j.k  - 


35xl0'6MeV  AZj 


■NiED-i-j'k(cos  ctg  )Aftk 


where 


Aftk  =  d<X>  sin  ddt}  =  27t(cosak+i  -  cos  aO 


Summing  over  incident  energies  and  pitch  angles  we  can  calculate  the  total  no.  of 
ion  pairs  cm-3  produced  the  jth  slab,  Qj,  i.e. 


5 


Qj  =  £  qij.k 
i.k 


Integrating  the  ion  pairs  cnr3  produced  over  altitude,  where  Laird’s  power  law 
spectrum  was  used  to  calculate  the  49  mono-energetic  fluxes,  gives  a  higher  ionization  for 
the  modified  Laird  method  by  about  a  factor  of  1.5  compared  to  the  results  of  the 
unmodified  Laird  method  for  day  293  of  1989.  (See  figure  4.  Here  the  circles  represent 
the  results  of  the  unmodified  Laird  method,  the  tiangles  represent  Jackman’s  results,  while 
the  squares  represent  the  results  of  the  modified  Laird  method.) 

Using  Jackman's  exponential  spectrum  to  calculate  the  49  mono-energetic  fluxes 
for  the  modified  Laird  method  gives  a  result  which  very  closely  resembles  Jackman’s. 
(See  figure  5.  The  circles,  triangles,  and  squares  represent  the  same  as  in  figure  4.) 
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#32  7-FEB-1991  16:19:19.35 

From:  6187 :: JACKMAN 

To:  KUPHSX: : ARMSTRONG, KUPHSX: : LAIRD 

CC: 

Subj:  Ion  Pair  Production  for  1989 

February  7,  1991 
Tom  and  Claude, 

I  thought  that  maybe  if  we  concentrated  on  just  one  day,  that  we  might  be 
able  to  track  down  our  differences.  The  day  I  have  chosen  to  give  you 
more  detailed  information  is  the  day  that  I  compute  maximum  ion  pair 
production  in  the  year  1989. 

Here  is  the  raw  data  that  you  sent  to  me  for  day  293  of  year  1989 
for  channels  Pi  -  Pll: 


1  IMP 8  1989  DAILY  AVG.  PROTON  FLUXES 

0YR=1989  DAY**  2  0  0  0  XSE=  34.65  YSE=  -5.36  ZSE=*  11.92 

YEAR  DAY  HH  MM  SS  P1AV  P2AV  P3AV  P4AV  P5AV 


,1989  293  0  0  0 

1 

0YR=1989  DAY=  2 
YEAR  DAY  HH  MM  SS 
,1989  293  0  0  0 


0.806E+04  0 . 755E+04  0.855E+04  0.559E+04 

IMP 8  1989  DAILY  AVG.  PROTON  FLUXES 
000  XSE=  34.65  YSE=  -5.36  ZSE=  11.92 

P7AV  P8AV  P9AV  P10AV 

273.  87.3  165.  52.2 


0 . 350E+04 


P11AV 

85.5 


I  derived  flux  values  {#  cm-2  s-1  sr-1  MeV-1)  from  this  data.  I  combined 
the  PI  and  P2  channels  and  also  the  P3  and  P4  channels.  The  energy  intervals, 
middle  energy  chosen,  channels  taken,  and  fluxes  computed  are  given  below 
for  day  293. 


Energy  intervals 
(MeV) 

Middle  energy  chosen 
(MeV) 

Channels  taken 

Flux  computed 

(#  cm-2  s-1  sr-1  MeV-1) 

0.29-0.96 

0.5 

PI  &  P2 

15,429. 

0.96-4.6 

2.0 

P3  &  P4 

2572.6 

4.6-15. 

9.8 

P5 

222.87 

15.-25. 

20. 

P7 

85.313 

25.-48. 

36.5 

P8 

11.861 

48.-96. 

72. 

P9 

10.742 

96.-145. 

120.5 

P10 

3.329 

190.-440. 

315. 

Pll 

1.069 

Here  is  some  information  on  fits  to  the  data  for  day  293. 

The  relationship  used  is  Flux  -  F  *  exp(-E/E0)  where  Flux  is 
in  units  of  #  cm-2  s-1  sr-1  MeV-1. 

The  values  of  F  and  E0  and  the  energy  bands  over  which  I 
have  applied  them  are  given  below: 


F  E0 

(MeV) 


Energy  band  applied 
(MeV) 


0.280E+05 
0 . 561E+03 
0 . 131E+02 
0.673E+01 


0 . 837E+00 
0 . 106E+02 
0 . 358E+03 
0 . 171E+03 


0.29-4.6 

4.6-25 

25-96 

96-440 


Next  I  give  you  output  from  my  energy  degradation  program.  I  give  height (km) , 
ion  pair  production  per  cm+3  per  sec,  and  ion  pair  production 
per  cm+3  per  hour. 

DAY-  293.00 

HEIGHT (KM)  Q(ION  PAIRS/ (CM+3  SEC))  Q(ION  PAIRS/ (CM+3  HOUR)) 


1.090E+02 

9 . 613E+01 

3 . 461E+05 

1 . 070E+02 

1 . 342E+02 

4 . 830E+05 

1 . 050E+02 

1 . 907E+02 

6 . 864E+05 

1.030E+02 

2 . 719E+02 

9.788E+05 

1 . 010E+02 

3.787E+02 

1 . 363E+06 

9 . 900E+01 

4 . 966E+02 

1.788E+06 

9.700E+01 

5 . 928E+02 

2 . 134E+06 

9 . 500E+01 

6 . 920E+02 

2 . 491E+06 

9 . 300E+01 

7 . 831E+02 

2 . 819E+06 

9 . 100E+01 

8 . 524E+02 

3 . 069E+06 

8 . 900E+01 

8 . 968E+02 

3 . 228E+06 

8.700E+01 

9 . 152E+02 

3.295E+06 

8.500E+01 

8 . 97  4E+02 

3.231E+06 

8 . 300E+01 

8 . 520E+02 

3 . 067E+06 

8 . 100E+01 

8 . 312E+02 

2 . 992E+06 

7 . 900E+01 

8 . 477E+02 

3 . 052E+06 

7 . 700E+01 

9.332E+02 

3.360E+06 

7 . 500E+01 

1.112E+03 

4 . 004E+06 

7 . 300E+01 

1 . 356E+03 

4 . 881E+06 

7 . 100E+01 

1 . 606E+03 

5.781E+06 

6. 900E+01 

1.766E+03 

6.359E+06 

6.700E+01 

1 . 907E+03 

6.865E+06 

6.500E+01 

2 . 018E+03 

7.266E+06 

6 . 300E+01 

2 . 076E+03 

7.473E+06 

6 . 100E+01 

2.070E+03 

7 . 453E+06 

5 . 900E+01 

2 . 047E+03 

7 . 370E+06 

5.700E+01 

1 . 960E+03 

7 . 055E+06 

5 . 500E+01 

1 . 877E+03 

€ . 758E+06 

5.300E+01 

1.820E+03 

6.554E+06 

5 . 100E+01 

1 . 966E+03 

7 . 078E+06 

4 . 900E+01 

2.303E+03 

8.290E+06 

4.700E+01 

2 . 684E+03 

9 . 662E+06 

4 . 500E+01 

3. 122E+03 

1 . 124E+07 

4 . 300E+01 

3 . 604E+03 

1.298E+07 

4 . 100E+01 

4 . 038E+03 

1 . 454E+07 

3.900E+01 

4 . 451E+03 

1 . 602E+07 

3 . 700E+01 

4 . 639E+03 

1 . 670E+07 

3 . 500E+01 

4 . 637E+03 

1 . 669E+07 

3 . 300E+01 

4.893E+03 

1.762E+07 

3 . 100E+01 

5 . 284E+03 

1 . 902E+07 

2 . 900E+01 

5 . 516E+03 

1.986E+07 

2.700E+01 

5 . 525E+03 

1 . 989E+07 

2.500E+01 

5 . 138E+03 

1 . 850E+07 

2 . 300E+01 

4 . 365E+03 

1 . 571E+07 

2 . 100E+01 

3 . 084E+03 

1 . 110E+07 

1 . 900E+01 

1 . 461E+03 

5.260E+06 

1 . 700E+01 

1 . 519E+02 

5 . 469E+05 

1 . 500E+01 

0 . 000E+00 

0 . 000E+00 

1 . 300E+01 
1 . 100E+01 
9 .  OOOE+OO 
7 . OOOE+OO 
5. OOOE+OO 
3. OOOE+OO 
1. OOOE+OO 
0.  OOOE+OO 


0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 


0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 


Here  is  the  standard  atmosphere  that  I  use: 


US.  STND. 
45  AVE 


ATMOS.  1976 


(km) 

Temp. 

Density (kg/m+3) 

110. 

240.0 

9 . 708E-08 

108. 

223.3 

1 . 381E-07 

106. 

212.9 

1 . 954E-07 

104. 

205.3 

2.769E-07 

102. 

199.5 

3.935E-07 

100. 

195.1 

5 . 604E-07 

98. 

191.7 

8 . 071E-07 

96. 

189.3 

1.162E-06 

94. 

187.7 

1 . 670E-06 

92. 

187.0 

2 . 393E-06 

90. 

186.9 

3 . 416E-06 

88. 

186.9 

4.875E-06 

86. 

186.9 

6.958E-06 

84. 

190.8 

9 . 694E-06 

82. 

194.7 

1 . 342E-05 

80. 

198.6 

1 . 846E-05 

78. 

202.5 

2 . 524E-05 

76. 

206.4 

3 . 431E-05 

74. 

210.4 

4 . 639E-05 

72. 

214.3 

6.237E-05 

70. 

219.6 

8.283E-05 

68. 

225.1 

1 . 092E-04 

66. 

230.5 

1.430E-04 

64. 

236.0 

1 . 861E-04 

62. 

241.5 

2 . 407E-04 

60. 

247.0 

3 . 097E-04 

58. 

252.5 

3 . 963E-04 

56. 

258.0 

5 . 045E-04 

54. 

263.5 

6 . 390E-04 

52. 

269.0 

8.056E-04 

50. 

270.7 

1 . 027E-03 

48. 

270.7 

1.317E-03 

46. 

266.9 

1.714E-03 

44. 

261.4 

2.259E-03 

42. 

255.9 

2 . 995E-03 

40. 

250.4 

3 . 996E-03 

38. 

244.8 

5.367E-03 

36. 

239.3 

7.258E-03 

34. 

233.7 

9 . 887E-03 

32. 

228.5 

1 . 356E-02 

30. 

226.5 

1.841E-02 

28. 

224.5 

2 . 508E-02 

26.  222.5  3 . 426E-02 

24.  220.6  4 . 694E-02 

22.  218.6  6.645L-02 

20.  216.7  8.891E-02 

18.  216.7  1.217E-01 

16.  216.7  1.665E-01 

14.  216.7  2.279E-01 

12.  216.7  3 . 119E-01 

10.  223.3  4 . 135E-01 

8.  236.2  5.258E-01 

6.  249.2  6 . 601E-01 

4.  262.2  8 . 194E-01 

2.  275.2  1 . 007E+00 

0.0  288.2  1.225E+00 

• 

Next  are  a  couple  of  simple  spectra  with  resultant  ion  pair  production. 

First,  an  exponential  fit:  Flux=  100  exp(-E/30)  where  E  is  in  MeV 

and  Flux  is  in  #  cm-2  s-1  sr-1  MeV-1. 

Results  from  my  energy  degradation  program  using  this  spectra: 


3.900E+01 

1 . 911E+03 

6.878E+06 

3.700E+01 

1 . 566E+03 

5 . 636E+06 

3. 500E+01 

1.210E+03 

4 . 357E+06 

3.300E+01 

B.300E+02 

2.988E+06 

3 . 100E+01 

5 . 134E+02 

1 . 848E+06 

2.900E+01 

2 . 783E+02 

1 . 002E+06 

2.700E+01 

1.232E+02 

4 . 436E+05 

2 . 500E+01 

4 . 878E+01 

1.756E+05 

2.300E+01 

1.371E+01 

4 . 936E+04 

2 . 100E+01 

3 . 189E+00 

1 . 148E+04 

1 . 900E+01 

3 . 363E-01 

1.211E+03 

1.700E+01 

0.000E+00 

0.000E+00 

1 . 500E+01 

O.OOOE+OO 

0.000E+00 

1 . 300E+01 

0.000E+00 

O.OOOE+OO 

1 . 100E+01 

0.000E+00 

O.OOOE+OO 

9 . 000E+00 

0.000E+00 

O.OOOE+OO 

7 . 000E+00 

0.000E+00 

O.OOOE+OO 

5 . 000E+00 

0.000E+00 

O.OOOE+OO 

3 . 000E+00 

O.OOOE+OO 

O.OOOE+OO 

1 . 000E+00 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

0.000E+00 

O.OOOE+OO 

Second,  a  power  law  fit:  Flux=  100,000.  *  (E/1) ** (-2)  where  E  is  in  MeV 
and  Flux  is  in  #  cm-2  s-1  sr-1  MeV-1. 

Results  from  my  energy  degradation  program  using  this  spectra: 

HEIGHT (KM)  Q(ION  PAIRS/ (CM+3  SEC))  Q(ION  PAIRS/ (CM+3  HOUR)) 


1 . 090E+02 

2 . 129E+03 

7 . 665E+06 

1 . 070E+02 

3 . 015E+03 

1 . 085E+07 

1 . 050E+02 

4 . 327E+03 

1.558E+07 

1 . 030E+02 

6 . 124E+03 

2.205E+07 

1 . 010E+02 

8 . 388E+03 

3 . 020E+07 

9 . 900E+01 

1 . 022E+04 

3 . 679E+07 

9.700E+01 

1 . 009E+04 

3 . 634E+07 

9 . 500E+01 

9.833E+03 

3 . 540E+07 

9 . 300E+01 

9 . 608E+03 

3 . 459E+07 

9 . 100E+01 

9.426E+03 

3 . 393E+07 

8 . 900E+01 

9 . 309E+03 

3 . 351E+07 

8 . 700E+01 

9 . 241E+03 

3 . 327E+07 

8 . 500E+01 

9 . 083E+03 

3 . 270E+07 

8 . 300E+01 

8 . 857E+03 

3.188E+07 

8 . 100E+01 

8 . 656E+03 

3.116E+07 

7 . 900E+01 

8.466E+03 

3 . 048E+07 

7 . 700E+01 

8 . 304E+03 

2 . 989E+07 

7 . 500E+01 

8 . 159E+03 

2 . 937E+07 

7 . 300E+01 

7 . 978E+03 

2.872E+07 

7 . 100E+01 

7 . 835E+03 

2.821E+07 

6 . 900E+01 

7 . 629E+03 

2.746E+07 

6.700E+01 

7 . 469E+03 

2 . 689E+07 

€ . 500E+01 

7 . 287E+03 

2.623E+07 

6 . 300E+01 

7 . 084E+03 

2.550E+07 

6 . 100E+01 

6 . 966E+03 

2 . 508E+07 

5 . 900E+01 

6 . 799E+03 

2.448E+07 

5.700E+01 

6 . 633E+03 

2.388E+07 

5 . 500E+01 

6 . 485E+03 

2 . 335E+07 

5 . 300E+01 

6 . 382E+03 

2.298E+07 

5 . 100E+01 

6.234E+03 

2.244E+07 

4 . 900E+01 

6.235E+03 

2 . 245E+07 

4.700E+01 

6.227E+03 

2.242E+07 

4 . 500E+01 

6 . 305E+03 

2.270E+07 

4 . 300E+01 

6 . 438E+03 

2 . 318E+07 

4 . 100E+01 

6 . 494E+03 

2 . 338E+07 

3 . 900E+01 

6 . 560E+03 

2 . 362E+07 

3 . 700E+01 

6 . 602E+03 

2 . 377E+07 

3 . 500E+01 

6 . 611E+03 

2 . 380E+07 

3 . 300E+01 

6 . 515E+03 

2 . 345E+07 

3 . 100E+01 

6 . 333E+03 

2 . 280E+07 

2 . 900E+01 

5 . 945E+03 

2 . 140E+07 

2.700E+01 

5 . 371E+03 

1 . 934E+07 

2 . 500E+01 

4 . 628E+03 

1 . 666E+07 

2 . 300E+01 

3 . 530E+03 

1 . 271E+07 

2 . 100E+01 

2 . 198E+03 

7 . 912E+06 

1 . 900E+01 

5 . 981E+02 

2 . 153E+06 

1 . 700E+01 

0. 000E+00 

0. 000E+00 

1 . 500E+01 

0. 000E+00 

0. 000E+00 

1 . 300E+01 

0. 000E+00 

0. 000E+00 

1 . 100E+01 

0. 000E+00 

0 .000E+00 

9 . 000E+00 

0. 000E+00 

0. 000E+00 

7 . OOOE+OO 

0. 000E+00 

0. 000E+00 

5 .  OOOE+OO 

0. 000E+00 

0. 000E+00 

3. 000E+00 

0. 000E+00 

0. 000E+00 

1. 000E+00 

0. 000E+00 

0. 000E+00 

0. 000E+00 

0. 000E+00 

0. OOOE+OO 

I  am  faxing  a  plot  of  the  proton  flux  versus  energy  for  day  293.  The 
data  are  represented  as  +'s  and  my  fit  as  a  solid  line.  The  discontinuities 
arise  when  the  e-folding  and  amplitude  values  change  for  the  different 
energy  bands . 

I  doubt  if  we  will  ever  compute  exactly  the  same  ion  pair  production,  but 
it  would  be  nice  to  be  within  a  factor  of  two  of  eachother. 


Thanks  again  for  spending  time  on  this. 

Sincerely, 

Charley  Jackman 


SPAN  address  -  PACF: : JACKMAN 
Phone  (301)286-8399 
FAX  (301)286-263 


Apendix  B 


Oct.  3,  1992 

TO:  C.  Jackman,  T.  P.  Armstrong,  C.  Laird,  T.  E.  Cravens,  G.  Dreschhoff,  E.  Zeller 
FROM:  Francis  Vitt 


After  having  resolved  the  differences  between  C  Jackman's  method  and  C.  Laird's 
and  T.  P.  Armstrong's  method  of  calculating  ionization  of  the  atmosphere,  to  bring  us  up 
to  date  on  the  computational  work  done  on  the  calculations,  the  following  is  an  outline  of 
the  Armstrong/Laird  method,  how  it  is  becoming  stream  lined,  and  some  proposed  goals  of 
the  problem. 

Claude  Laird's  computational  procedure  is  summarized  in  the  following  diagram 
(figure  1).  Here  we  are  assuming  that  the  energy  spectrum  of  the  proton  fluxes  can  be 

sg£  =  A*E-5  dR 

validly  described  by  the  power  law  dE  ,  where  dE  is  the  proton  differential  flux 

(number  of  protons  per  square  cm  sr  sec  MeV)  and  8  is  the  spectral  exponent 


The  differential  flux  for  each  energy,  E,  (which  are  10  logarithmic  midpoints  of 
IMP8  proton  channel  energies,  contained  in  the  data  file  DA  TUM DAT),  and  the  spectral 
exponent, 8 ,  for  each  of  the  9  energy  ranges  (contained  in  the  data  file  ERANGEDAT) 
between  the  above  10  logarithmic  spaced  energies,  are  calculated  using  the  program 
NTOTAL1  FOR  which  reads  data  from  a  daily  average  IMP8  tape  and  outputs  the  above 
calculated  quantities  to  the  tapes  called  TOTAL  TAPES. 


Reading  from  the  TOTAL  TAPES ,  the  program  NTOTREC1FOR,  using  the 
subroutine  NP ARTICLES. FOR,  performs  the  integration 

m.  -  24  hrs  *3.6QQscc  *  I 
1  J—  J  dis 


day 


hr 


,  where  Ej  (i  =  1,2,. ..,49)  are  logarithmic  spaced 
energies  ranging  from  0.38  MeV  to  289  MeV,  giving  49  mono-energetic  fluxes  (number  of 
protons  per  square  cm  sr  sec)  incident  on  the  upper  atmosphere  with  incident  energies  E?  = 
YEj*Ej+i, , which  are  output  in  the  data  file  49DAILYDAT.  NTOTREC1  FOR  also 
outputs  the  data  file  TOTALDAT.  which  contains  the  differential  flux  and  exponent  data 
for  the  period  considered. 


The  program  NPAMYMUTDAYFOR,  which  inputs  the  49  mono-energetic  fluxes 
contained  in  49DAJLYDAT,  follows  these  fluxes  through  the  atmosphere  as  energy 


deposited  in  each  of  the  5  km  slabs  is  calculated  using  a  table  look  up  method  using  range 
and  energy  data  contained  the  data  files  RANGE RNG  and  RANGEFNG.  Subroutine 
SLAB. FOR  returns  to  the  program  NPAMY MUTDAY BOR  the  equivalent  thicknesses  of 
the  slabs  at  sea  level  using  monthly  average  atmospheric  densities  at  the  top  and  bottom  of 
each  atmospheric  slab,  contained  in  the  data  file  ATMOS2DAT.  NPAMYMUTDAY  FOR 
outputs  the  total  energy  deposited  per  cubic  cm  in  each  slab  and  total  number  of  ion  pairs 
per  cubic  cm  produced  in  each  slab  produced  by  the  protons  incident  on  the  upper 
atmosphere  for  each  day  in  the  data  file  NDELTAEDAYDAT. 

A  similar  process  is  used  to  calculate  the  total  energy  deposited  and  the  ion  pairs 
produced  in  the  5  km  slabs  by  the  alpha  particles.  Summing  the  two  gives  the  total  number 
of  ion  pairs  per  cubic  cm  produced  by  both  the  alpha  particles  and  protons  in  each  slab  each 
day. 

Figure  1: 


The  above  process  has  been  stream  lined  to  calculate,  using  input  data  from  an  IMP 
8  daily  average  tape,  the  total  number  of  ion  pairs  produced  by  both  alpha  particles  and 
protons  per  cubic  cm  in  each  5  km  slab  for  each  day,  which  is  output  in  the  data  file 
IONPAIRSnnDAT.  In  short,  to  summarize  the  stream  lined  process: 


Input: 


IMP  8  daily  average  tape. 

Starting  and  end  dates  of  the  period  desired. 

Output: 

IONPAIRSnnDAT  (nn  equals  the  year  minus  1900,  e.g.  for  1990  nn  =  90) 

The  following  diagram  (figure*.  2)  summarizes  how  the  process  has  been  stream 
lined  • 

Program  NTOTAL1  FOR  has  be  ,n  converted  into  the  subroutine 
NTOTSUBFOR,  which  is  driven  by  the  program  IONJPRSFOR.  The  user  inputs  into 
IONPRSFOR  the  starting  and  ending  year  and  day  of  year  of  the  dates 
IONPAIRSnnDAT  is  wanted. 

ION  PRSFOR  calls  subroutine  NTOT_SUBFOR  once  each  year  of  the  desired 
dates.  Subroutine  NTOTSUBFOR  reads  data  from  the  IMP8  tape,  calls  subroutines 
NPARTS UBFOR,  NPARTAJ  SUB  FOR,  NPAMSUBFOR,  NPAMALJUBFOR, 
and  outputs  data  to  IONPAIRSnnDAT  once  each  day  of  each  year. 

Subroutine  NPART  S UB FOR  (subroutine  NPARTALSUB FOR)  is  the 
modified  above  subroutine  N PARTICLES  FOR,  which  calculates  the  49  mono-encrgetic 
proton  fluxes  (35  mono-energetic  alpha  particle  fluxes). 

Subroutine  NPAM  SUBFOR  (subroutine  NPAMALJUBFOR), which  is  the 
modified  program  NPAMYMUTDAYFOR,  calculates  the  total  energy  deposited  and 
number  of  ion  pairs  per  cubic  cm  produced  in  each  slab  by  the  49  (35)  mono-energetic 
fluxes  by  using  the  subroutine  DEPOSIT  FOR  {DEPOSIT AL.FOR). 

Subroutine  DEPOSITFOR  (subroutine  DEPOSIT ALFOR)  integrates  over  pitch 
angle  from  0  to  pi/2  (assuming  an  isotropic  distribution  of  pitch  angle,  and  symmetric 
zenith  angle)  to  calculate  energy  deposited  by  the  table  look  up  method  using  data  files 
RANGE.ENG  and  RANGE.RNG  {AIR_ALPHAS.ENG  and  AJR_ALPHASRNG)  by 
each  of  the  mono-energetic  fluxes. 

Subroutine  SLAB2JUBFOR,  modified  version  of  the  above  subroutine 
SLAB2FOR,  returns  equivalent  sea  level  thickness  to  the  subroutines  NPAMJUBFOR 


and  NPAMALJS UB FOR  for  each  date  using  monthly  average  densities  at  the  top  and 
bottom  of  each  slab,  which  is  contained  in  the  data  file  MS _SLAB SnnJDAT  (nn  =  year  - 
1900). 


The  data  files  MS_SlABSnnDAT  were  created  using  MSIS  (Mass  Spectrometer 
and  Incoherent  Scatter  model)  and  MSISE  (the  extended  MSIS  model)  atmospheric 
models. 

The  MSIS  model  has  a  lower  limit  of  80  km,  the  MSISE  model  is  for  altitudes 
greater  than  0  km  (the  latter  is  the  extended  version  of  the  former,  but  is  not  recommended 
over  the  former  for  altitudes  greater  than  85  km). 

These  models,  which  were  obtained  from  NSSDC  (NASA's  National  Space 
Science  Center),  require  the  parameters  day  of  year,  universal  time,  altitude,  geodetic 
latitude  and  longitude,  local  apparent  solar  time,  three  month  average  of  10.7cm  solar  radio 
flux,  daily  10.7cm  solar  radio  flux  for  previous  day,  and  daily  AP  magnetic  activity  index. 
The  daily  10.7cm  solar  radio  flux  and  the  daily  AP  magnetic  activity  index  were  obtained 
from  Solar  Geophysical  Data. 


Figure  2: 


The  results  of  running  the  process  using  the  atmospheric  model  which  Claude  Lai: 1 
has  used  has  been  compared  to  the  results  of  using  the  MSIS  and  MSISE  atmospheric 
models.  The  largest  discrepancies  appear  to  occur  during  periods  of  relatively  large  soiai 


activity.  To  illustrate  this,  figures  3  and  4  are  contour  plots  of  the  log  of  number  of  ion 
pairs  per  cubic  cm  per  hr  by  protons  and  alpha  particles  from  the  results  of  running  the 
process  using  Claude  Laird's  atmospheric  model  for  years  1990  and  1991,  while  figures  5 
and  6  are  the  result  of  using  the  MSIS  models.  Figures  7  and  8  are  contour  plots  of  the  log 
of  absolute  value  of  the  difference  in  the  results  of  using  the  two  models  for  199C  .J 
1991.  As  you  can  see,  the  larger  discrepancy  occurs  during  the  periods  of  when  solar 
events  have  occurred.  This  is  probably  a  result  of  the  MSIS  models  use  the  daily  10.7cm 
solar  radio  flux  and  the  daily  AP  magnetic  activity  index,  which  are  high  during  these 
periods. 

A  suggested  goal  of  the  calculation  is  to  link  the  input  ion  pair  distribution  with 
vertical  transport  of  the  ions  and  their  life  times.  That  is,  after  calculating  the  total  number 
of  ion  pairs  per  cubic  cm  produced  in  each  slab  for  a  given  day,  we  need  to  add  to  the  total 
the  number  of  ion  pairs  per  cubic  cm  which  do  not  recombine,  and  are  transported  into  each 
slab  from  slabs  above  and  below  from  previous  days.  So  that  we  are  calculating  the  total 
number  ion  pairs  per  cubic  cm  in  each  slab  each  day,  not  just  the  ionization  produced  each 
day. 


ABS(DIFF  ION  PAIRS/CUBIC  CM-HR) 
(PROTONS  +  ALPHAS) 


(PROTONS  +  ALPHAS) 


August  6,  1992 


Regarding  the  calculation  of  ion  pairs  produced  by  energetic  protons  incident  on  the 
upper  atmosphere  using  IMP8  data  which  has  been  done  by  Claud  Laird  and  Charlie 
Jackman,  daily  average  IMP8  data  for  day  293  of  1989  has  been  used.  Integration  of  ion 
pairs  per  cubic  cm  over  the  altitude  shows  that  Jackman's  method  tends  to  give  more 
ionization  than  Laird's  method,  by  about  a  factor  of  two.  (See  figure  1.  Here  the  circles 
represent  Laird's  results,  triangles  represent  Jackman's  results.)  To  find  the  discrepancy 
between  the  two  methods  the  following  comparisons  were  done.  Modifications  were  done 
to  the  Laird  method  to  attempt  to  obtain  results  that  match  Jackman's. 

Comparison  of  the  two  methods: 


*  Both  assume  that  the  number  of  ion  pairs  produced  equals  the  energy  deposited 
by  the  protons  in  the  atmosphere  divided  by  35  eV. 

^-  =  A*E-5  dK 

*  Laird  uses  the  power  law  energy  spectrum  ,  where  JE  is  the 

differential  flux  (no.  of  particles  *cnr2*sr1*sec.*1*  MeV'1),  A  is  the  flux  constant,  5  is 
the  spectral  exponent  A  and  5  are  determined  over  nine  energy  ranges  to  fit  the  IMP8  data. 
This  gives  a  piece  wise  continuous  spectrum.  (See  figure  2.) 


^  =  F*exp(-E-) 

*  Jackman  uses  an  exponential  energy  spectrum  given  by  dE  Eo  ,  where 

F  and  Eo  are  determined  over  only  four  energy  ranges  to  fit  the  the  IMP8  data.  This  gives 
relatively  large  discontinuities  in  the  energy  spectrum.  (See  figure  3.) 


*  Laird's  energy  spectrum  ranges  from  0.38  MeV  to  289  MeV. 


*  Jackman’s  energy  spectrum  ranges  from  0.29  MeV  to  440  MeV. 


*  Laird's  method: 


Integration  of  the  energy  spectrum  gives 

N  _  24_hrs  *3600sec  *  I  dN^g 
1  day  hr  I  dE 


=  total  no.  of  particles  cnr2  sr 1  incident  on  the  upper 
atmosphere  for  the  day, 

where  Ej  (i  =  1,2,. ..,50)  are  logarithmically  spaced  energies  ranging  from  0.38  MeV  to 
289  MeV.  This  gives  49  mono-energetic  fluxes  with  incident  energy  E?  =  VEi*Ei+K 


These  fluxes  are  followed  vertically  through  the  atmosphere,  which  is  divided  into 
5  km  slabs,  as  energy  deposited  in  each  slab  is  calculated,  using  a  table  look  up  method. 

Hence  calculating  energy  deposited  cm-2  sr1.  Then  Laird  multiplies  by 
( - 1 - 1 - 2JLS) 

35xl0'6MeV  5x1 05^  3  to  get  the  units  of  no.  of  ion  pairs  cnr3. 


*  Jackman's  method: 


Jackman  assumes  an  isotropic  distribution  of  pitch  angle  and  azimuthal  symmetry. 
Jackman  calculates  no.  of  ion  pairs*cnr3*sec.'1  using  the  expression 

qM-Pttf  fg^dEdn 


m 

where  dE  is  the  differential  flux  (no.*cnr2*sr1*sec.'1*MeV'1),  w  =  35  eV  per  ion  pair, 

dE 

is  the  atmospheric  density  at  the  altitude  z,  dx  is  the  energy  deposited  per  gm. 


Let  Ed  =  dE  be  the  energy  deposited  by  a  particle  of  energy  E  with  a  pitch  angle  a 
in  the  segment  zt>,  3D  =  width  of  the  atmospheric  layer  in  gm*cnr2  centered  on  the  altitude 
z. 


d  (Eq  gc  z) 

If  dE  ’  ’  =  flux  of  particles  with  incident  energy  Eo  at  a  pitch  angle  a  at  the 
altitude  z  is  known,  then  9i,j(^>z)  =  ionization  rate  by  a  flux  of  particles  in  the  energy 


Eo  -  4P-  to  Eo  + 

interval  2  * 


a_4fittoa  +  £& 

2  ,  and  in  the  pitch  angle  interval  2  2  over  the 


atmospheric  distance  zd  centered  on  the  altitude  z  can  be  calculated  with  the  expression 


Ja.AS 

2 


sjg<Eo,a,z) 

Ed^E - 

w  zd  seca 


qitj(A,z)  =2n 


sina  da  dE 


dii(E0)a,z) 

Assuming  Ed  and  dE v  ’  ’  ’  change  slowly  enough  within  the  pitch  angle  interval 


we  can  write 


^<Eo,a,z)  Ed  dE  I  cos  a  sin  a  da 

Z).i£  Ja.i 

/*>*¥ 

=  ^  [sin2(a+^)  -  sin2(a*  *&■)]  ^{Eo.a.z)  ed  dE 

/fc-AE 


To  calculate  the  energy  deposited  Jackman  uses  the  expression  for  the  range  as 
R(E)=A‘EB,thusED  =  E-(-a!fS  +  EB)i 


Hence,  the  ionization  rate  at  the  altitude  z  is 
q(z)  =  p(z)  X  qi.j(A,z) 


=  no.  of  ion  pairs  cm*3  sec.*1 


Modifications  to  the  Laird  method: 


Assume  an  isotropic  distribution  of  pitch  angle  between  0  and  2  and  azimuthal 

2 L 

symmetry.  Let  ak  be  50  evenly  spaced  angles  between  0  and  2  ,  the  fluxes  with  pitch 
angle  between  and  ak+i  have  the  pitch  angle  ak  such  that cos  ak  =  Vcos  ajc*cos  ak+i 

Let  qi-j.k  =  the  no.  of  ion  pairs*cm*3  produced  by  the  flux  with  incident  energy  E? 
in  the  J*  slab  centered  at  altitude  zj  with  pitch  angle  ak.  So, 


qi.j,k  = 


35xlO*6MeV 


f 

L  AXj>k  # 


where  ED,i.j,k  is  the  energy  deposited  by  a  proton  with  incident  energy  E?,  pitch  angle  ak, 
in  the  1th  slab,  is  the  effective  thickness  of  the  1th  slab,  i.e.  =  AZj*secakj  AZj 
is  the  thickness  of  the  slab  (5km).  Nj  (for  i  =  1,2, ...49)  is  the  mono-energetic  flux  with 
incident  energy  E?.  And  AQ  is  the  solid  angle  between  ak  and  ak+i  as  0  varies  from  0  to 


Assuming  Njt  ED.i.jJc,  and  do  not  change  significantly  as  pitch  angle  varies 
between  «ic  and  <*k+i,  we  can  write 


9i,j.k  — 


_L 


SSxlO^MeV  AZj 


N,- — ‘3'k-(cos  a°  )AQic 


where 


AQjc  =  j  dO  I  sin  ddd  =  27t(cosak+i  -  cos  a*) 

J* k 


Summing  over  incident  energies  and  pitch  angles  we  can  calculated  the  total  no.  of 
ion  pairs  cm-3  produced  the  the  1th  slab,  Qi,  i.e. 

Qj  =  X  9i.j.k 

i.k 


Integrating  the  ion  pairs  cm'3  produced  over  altitude,  where  Laird's  power  law 
spectrum  was  used  to  calculate  the  49  mono-energetic  fluxes,  gives  gives  a  higher 
ionization  for  the  modified  Laird  method  by  about  a  factor  of  1.5  compared  to  the  resuits  of 
the  unmodified  Laird  method  for  day  293  of  1989.  (See  figure  4.  Here  the  circles 
represent  the  results  of  the  unmodified  Laird  method,  the  triangles  represent  Jackman's 
results,  while  the  squares  represent  the  results  of  the  modified  Laird  method.) 

Using  Jackman’s  exponential  spectrum  to  calculate  the  49  mono-energetic  fluxes 
for  the  modified  Laird  method  gives  a  result  which  very  closely  resembles  Jackmam’s. 

(See  figure  5.  The  circles,  triangles,  and  squares  represent  the  same  as  in  figure  4.) 
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Apendix  C 


June  8,  1993 
MEMO 

FROM;  Francis  Vitt 


During  solar  proton  events  large  number  of  energetic  protons,  which  are  detected 
by  IMP8,  are  released  by  the  sun.  By  using  an  energy  deposition  program  we  are  able  to 
calculate  ion  pair  production  rates  versus  altitude  from  IMP8  daily  averaged  proton  count 
rates.  The  energy  deposited  is  calculated  as  the  protons  are  allowed  to  propagate  down 
through  the  extended  MSIS  model  atmosphere  (ref.  Hedin  A.E.,  Extension  of  the  MSIS 
Thermosphere  Model  into  the  Middle  and  Lower  Atmosphere,  J.G.R.,  96,  A2,  1 159, 
1991). 


These  ion  pair  production  rate  profiles  have  been  compared  to  the  ion  pair 
production  rates  calculated  by  Nicolet  in  1975  due  to  galactic  cosmic  rays,  i.e.  protons 
from  outside  the  solar  system,  (ref.  Nicolet,  M.,  On  the  production  of  Nitric  Oxides  by 
Cosmic  Rays  in  the  Mesosphere  and  Stratosphere ,  Planet  Space  Sri.,  23,  637, 1975).  It 
was  found  that  the  ion  pair  production  rates  calculated  using  the  IMP8  data  fell  short  of  the 
ion  pair  production  rates  calculated  by  Nicolet  by  about  an  order  of  magnitude  in  the  lower 
altitudes.  This  may  be  caused  by  the  fact  that  IMP8  can  not  detect  protons  of  energies 
greater  than  440  MeV  and  the  energies  of  the  galactic  ray  protons  are  much  greater  than 
this.  This  is  remedied  by  adjusting  the  IMP8  data  by  subtracting  the  back  ground  from  the 
count  rates,  using  the  adjusted  data  to  the  calculate  ion  pair  production  rates  with  the  energy 
deposition  program,  then  adding  to  these  ion  pair  production  rates  the  rates  calculated  by 
Nicolet,  i.e.  the  ion  pair  production  rates  due  to  the  galactic  cosmic  rays  are  added  to  the 
adjusted  production  rates  due  to  the  solar  proton  events. 

With  an  estimated  ratio,  obtained  from  the  literature,  of  the  ion  pairs  produced  to 
odd  nitrogen  species  produced  gives  a  production  rate  of  odd  nitrogen.  Using  the 
production  rates  due  to  solar  proton  events  as  odd  nitrogen  production  rate  coupled  with 
chemical  production  rates  and  loss  rates,  coupled  with  transport  processes  will  give  a  time 
dependent  model  of  the  odd  nitrogen  profiles.  Currently,  a  one  dimensional  diffusion 
problem  with  appropriate  source  and  loss  terms  in  the  continuity  equation  is  being  solved. 


